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C ■ ABSTRACT 

: 

We find the distribution of coalescence times, birthrates, spatial velocities, and 
lO ' subsequent radial offsets of coalescing neutron stars (NSs) in various galactic potentials 

accounting for large asymmetric kicks introduced during a supernovae. The birthrates 
of bound NS-NS binaries are quite sensitive to the magnitude of the kick velocities but 
are, nevertheless, similar (~ 10 per Galaxy per Myr) to previous population synthesis 
studies. The distribution of merger times since zero-age main sequence is, however, 
' relatively insensitive to the choice of kick velocities. With a median merger time of 

I ~ 10® yr, we find that compact binaries should closely trace the star formation rate 

i in the Universe. 

' In a range of plausible galactic potentials (with Afgaiaxy <; 3 x 10^°Mq) the median 

OO . radial offset of a NS-NS mergers is less than 10 kpc. At a redshift of z = 1 (with 

0\ • = 65 km s-i Mpc-i and n = 0.2), this means that half the coalescences should 

fH ^ occur within '--^ 1.3 arcsec from the host galaxy. In all but the most shallow potentials, 

ninety percent of NS-NS binaries merge within 30 kpc of the host. We find that 
although the spatial distribution of coalescing neutron star binaries is consistent with 
the close spatial association of known optical afterglows of gamma-ray bursts (GRBs) 
with faint galaxies, a non-negligible fraction (^ 15 percent) of GRBs should occur well 
outside (^ 30 kpc) dwarf galaxy hosts. Extinction due to dust in the host, projection 
of offsets, and a range in interstellar medium densities confound the true distribution 
of NS-NS mergers around galaxies with an observable set of optical transients/galaxy 
offsets. 
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1 INTRODUCTION 

The discovery of an X-ray afterglow by BeppoSAX (Costa et 
al. 1997; Piro, Scarsi, & Butler 1995) and subsequently an 
optical transient associated with gamma-ray burst (GRB) 
970228 (van Paradijs et al. 1997) led to the confirmation 
of the cosmological nature of GRBs (Metzger et al. 1997). 
The broadband optical afterglow has been modeled rela- 
tively successfully (Meszaros & Rees 1993; Wijers, Rees, & 
Meszaros 1997; Waxman 1997; Waxman, Kulkarni & Frail 
1998) as consistent with an expanding relativistic fireball 
(Rees & Meszaros 1994; Paczyhski & Rhoads 1993; Katz 
1994; Meszaros & Rees 1997; Vietri 1997; Sari, Piran, & 
Narayan 1998; Rees & Meszaros 1998). 

Still, very little is known about the nature of the progen- 
itors of GRBs, and, for that matter, their hosts. Broad-band 



fluence measures and the known redshifts of some bursts 
implies a minimum (isotropic) energy budget for GRBs of 
~ 1052-53 gj.gg (Metzger et al. 1997; Kulkarni et al. 1998; 
Djgorovski et al. 1998). The log A''-log P brightness distribu- 
tion, the observed rate, A'', of bursts above some flux, P, ver- 
sus flux, indicates a paucity of dim events from that expected 
in a homogeneous, Euclidean space. With assumptions of 
a cosmology, source evolution and degree of anisotropy of 
emission, the log A'^-log P has been modeled to find a global 
bursting rate. Assuming the bursts are non-evolving stan- 
dard candles Fenimore &c Bloom 1995 found ~ 1 bursts 
events per galaxy per Myr (GEM) to be consistent with 
the observed log A^-Iog P. More recently, Wijers et al. 1998 
(see also, Totani et al. 1998; Lipunov, Postnov, & Prokhorov 
1997) found the same data consistent with GRBs as stan- 
dard candles assuming the bursting rate traces the star- 
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formation rate (SFR) in the Universe; such a distribution 
implies a local burst rate of ~ 0.001 GEM and a stan- 
dard peak luminosity of Lo = 8.3 x 10^^ erg (Wijers 
et al. 1998). 

Given the energetics, burst rate and implied fluences, 
the coalescence, or merger, of two bound neutron stars (NSs) 
is the leading mechanism whereby gamma-ray bursts are 
thought to arise (Paczyiiski 1986; Goodman 1986; Eichler 
et al. 1989; Narayan, Paczyriski, & Piran 1992). One quan- 
tifiable prediction of the NS-NS merger hypothesis is the 
spatial distribution of GRBs (and GRB afterglow) with re- 
spect to their host galaxies. Conventional wisdom, using the 
relatively long lived Hulsc- Taylor binary pulsar as a model, 
is that such mergers can occur quite far (> 100 kpc) out- 
side of a host galaxy. Observed pulsar (PSll) binaries with 
a NS companion provide the only direct constraints on such 
a populations, but the observations are biased both towards 
long lived systems, and systems that are close to the Galac- 
tic plane. 

The merger rate of NS-NS binaries has been discussed 
both in the context of gravitational wave-detection and 
GRBs (eg. Phinney 1991; Narayan et al. 1991; Tutukov 
& Yungelson 1994; Lipunov et al. 1995). Recently Fryer, 
Burrows, and Benz (1997), Lipunov, Postnov, & Prokhorov 
(1997), Portegies Zwart & Spreeuw (1996) studied the ef- 
fect of asymmetric kicks on birthrates of NS-NS binaries, 
but did not quantify the spatial distribution of such bina- 
ries around their host galaxies. Tutukov & Yungelson (1994) 
discussed the spatial distribution of NS-NS mergers but ne- 
glected asymmetric kicks and the effect of a galactic po- 
tential in their simulations. Only Zwart & Yungelson (1998) 
have discussed the maximum travel distance of merging neu- 
tron stars including asymmetric supernovae kicks. 

It is certainly of interest to find the rate of NS-NS co- 
alescences ab initio from population synthesis of a stellar 
population. This provides an estimate of beaming of GRBs, 
assuming they are due to NS-NS mergers, and hence an es- 
timate of probable frequency of gravitational wave sources, 
providing a complementary rate estimate to those of Phin- 
ney 1991 and Narayan et al. 1991, which are based on long 
lived NS-PSR pairs only and are very conservative. It also 
provides an estimate of how the NS NS mergers trace the 
cosmological stax formation rate (SFR) of the Universe, if 
mean formation rates and binarity of high mass stars axe 
independent of star formation environments such as metal- 
licity. 

Here we concentrate on estimating the spatial distribu- 
tion of coalescing NS-NS binaries around galaxies. To do 
so, both the system velocity and the interval between for- 
mation of the neutron star binary and the merger through 
gravitational radiation is found by simulation of binary sys- 
tems in which two supernovae occur. We explore the effects 
of different asymmetric kick amplitudes, and the resultant 
birthrates and spatial distribution of coalescing NS-NS bi- 
naries born in different galactic potentials. 

In section 2 we briefly outline the prescription for our 
Monte Carlo code to simulate bound binary pairs from an 
initial population of binaries by including the effect of asym- 
metric supernovae kicks. In section 3 we outline the integra- 
tion method of NS-NS pairs in various galactic potentials. 
Section 4 highlights the birthrates and spatial distributions 
inferred from the simulations. Section 5 concludes by discus- 



sion the implications and predictions for gamma-ray burst 
studies. 



2 NEUTRON STAR BINARY POPULATION 
SYNTHESIS 

We used a modified version of the code created for binary 
evolution by Pols (Pols & Marinus 1994) taking into ac- 
count the evolution of eccentricity through tidal interaction 
and mass transfer before the first and second supernova, and 
allowing for an asymmetric kick to both the NS during su- 
pernova. The reader is referred to Pols & Marinus (1994) for 
a more detailed discussion account of the binary evolution 
code. 

2.1 Initial Conditions and Binary Evolution 

In general, the evolution of a binary is determined by the 
initial masses of the two stars (mi, rn-z), the initial semi- 
major axis (oo) and the initial eccentricity (eo) of the binary 
at zero-age main sequence (ZAMS). We construct Monte 
Carlo ensembles of high-mass protobinary systems (with pri- 
mary masses between 4Mq and lOOM©) by drawing from 
an initial distribution of each of the four parameters as pre- 
scribed and motivated in Portegies Zwart & Verbunt (1996). 
We treat mass transfer and common-envelope (CE) phases 
of evolution as in Pols & Marirms (1994). CE evolution is 
treated as a spiral-in process; we use a value of a = 1 for 
the efficiency parameter of conversion of orbital energy into 
envelope potential energy, see eq. [17] of Pols & Marinus 
(1994). We treat circularization of an initially eccentric or- 
bit as in Portegies Zwart & Verbunt (1996). 

During detached phases of evolution we assume that 
mass accreted by the companion is negligible so that aMtot 
— constant. Mass lost by the binary system in each suc- 
cessive time step results in a change in eccentricity accord- 
ing to the sudden mass loss equations (see, for example, 
eqns. [A.21] and [A.24] of Wettig & Brown 1996). We ignore 
the effect of gravitational radiation and magnetic braking in 
the early stages of binary evolution. 

The simple approximation of the 4-parameter distribu- 
tion function, albeit rather ad hoc, appears to adequately 
reproduce the observed population of lower mass stars in 
clusters (eg. Pols & Marinus 1994). The effect on the dis- 
tribution of NS NS binaries after the second supernova by 
variation of the 4-paramcter space is certainty of interest, 
but we have used the canonical values. A fair level of ro- 
bustness is noted in that varying the limits of the initial 
distributions of Oo and Co does not the change the implied 
birthrates of bound NS binaries nearly as much as plau- 
sible variation in asymmetric kick distribution. This effect 
was noted in Portegies Zwart & Spreeuw 1995 and Portegies 
Zwart & Verbunt 1995. 

2.2 Asymmetric Supernovae Kicks 

Several authors (eg. Paczyhski 1990; Narayan & Ostriker 
1990; Lyne & Lorimer 1994; Cordes & Clicrnoff 1997) have 
sought to constrain the distribution of an asymmetric kick 
velocities from observations of isolated pulsars which are the 
presumed by-products of type II supernovae. Even careful 
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modeling of the selection effects in observing such pulsars 
has yielded derived mean velocities that differ by nearly an 
order of magnitude. It is important here to use a good es- 
timate for the actual physical impulse (the "kick velocity") 
the neutron stars receive on formation. The observed distri- 
bution of pulsar velocities does not reflect the kick distribu- 
tion directly as it includes the Blaauw kick (Blaauw 1961) 
from those pulsars formed in binaries, and selection effects 
on observing both the high and low speed tail of the pulsar 
population (eg. Hartman 1997). Hansen & Phirmcy (1997) 
found that the observed distribution is adequately fit by a 
Maxwellian velocity distribution with Ukick = 190 km s~^ 
(which corresponds to a 3-D moan velocity of 300 km s~^). 
Since it is not clear that pulsar observations require a more 
complicated kick-velocity distribution, we chose to adopt a 
Maxwellian but vary the value of crkick. 

When a member of the binary undergoes a supernovae 
we assume the resulting NS receives a velocity kick, Vk, 
drawn from this distribution. Although the direction of this 
kick may be coupled to orientation of the binary plane, we 
choose a kick with a random spatial direction, since there is 
no known correlation between the kick direction or magni- 
tude and the binary parameters. 

With an angle a between the velocity kick and the rel- 
ative velocity vector, v, of the two stars. Then, following 
earlier formulae (e.g. Portegies Zwart & Verbunt 1996; Wet- 
tig & Brown 1996), the new-semi major axis of the binary 
is. 



'f'^ + i^k + 2t;wfe cos a 
G(Mns + M2) 



(1) 



where r is the instantaneous distance between the two stars 
before SN, M2 is the mass of the companion (which may 
already be a NS), and Mns = IAMq is the mass of the 
resulting neutron star. We neglect the effects of supernova- 
shell accretion on the mass of the companion star. If a' is 
positive, the new eccentricity is 



a'G(MNs + M2) 



1/2 



(2) 



where the resultant relative velocity is Vr = v+Vk- Assuming 
the kick directions between successive SN are independent, 
the resulting kick to the bound system (whose magnitude 
is given by equation [2.10] of Brant & Podsiadlowski 1995) 
is added in quadrature to the initial system velocity to give 
the system velocity (vsys)- 

To produce 1082 bound NS-NS binaries with a Hansen 
& Phinncy kick velocity distribution and initial conditions 
described above, we follow the evolution of 9.7 million main 
sequence binaries which produce a total of ~ 1 million neu- 
tron stars through supernovae. Assuming a supernova rate 
of 1 per 40 years (Tammann et al. 1994) and 40% binary 
fraction (as in Portegies Zwart & Spreeuw 1996), we find 
an implied birthrate of NS-NS binaries by computing the 
number of binaries with SN type II per year and multiply 
by the ratio of bound NS-NS systems to SN type II as found 
in the simulations. We neglect the (presumed small) contri- 
bution of other formation channels (eg. three-body interac- 
tions) to the overall birthrate of NS-NS binaries. The im- 
plied birthrate of NS-NS binaries from various kick-velocity 
magnitudes are given in table 2. 
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Figure 1. The distribution of orbital parameters (period and 
eccentricity) after the second supernovae for bound NS-NS pairs. 
Prom left to right, arc lines of constant merger time after second 
SN (10'', 10**, 10^" yrs). Tlie parameters of the observed NS pairs 
1913-1-16 (Taylor & Weisberg 1989), 1534-1-12 (Wolszczan 1991), 
and 2303-1-46 (Taylor & Dewey 1988) arc marked with triangles. 
With an observational bias towards long-lived systems, clearly 
the observed PSR-NS systems are not indicative of the true NS 
binary distribution (see section 4.1). 



3 EVOLUTION OF BINARIES SYSTEMS IN A 
GALACTIC POTENTIAL 

The large-scale dynamics of stellar objects are dominated 
by the halo gravitational potential while the initial distribu- 
tion of stellar objects is characterized by a disk scale length. 
We take the disk scale and halo scale to vary independently 
in our galactic models. We assume the NS-NS binaries axe 
born in an exponential stellar disk, with birthplace drawn 
from randomly from mass distribution of the disk. The ini- 
tial velocity is the local circular velocity (characterized but 
the halo) plus Vsys added with random orientation. 

We then integrate the motion of the binary in the galac- 
tic potential assuming a Hernquist (1990) halo; we ignore the 
contribution of the disk to the potential. We assume scale 
lengths for the disk and halo, the disk scale (raisk) deter- 
mines the disk distribution, the halo scale length (rbreak) and 
circular velocity (wdrc) determine the halo mass (see table 
1). The movement of the NS-NS binaries on long time-scales 
is sensitive primarily to the depth of the galactic potential 
(here assumed to be halo dominated) and how quickly it 
falls off at large radii. Assuming isothermal halos instead 
of Hernquist profiles would decrease the fraction of NS-NS 
pairs that move to large galactocentric radii, but the dif- 
ferences in distribution are dominated by the true depth of 
the halo potentials in which the stars form rather than their 
density profiles at large radii. 

We use a symplectic leapfrog integrator to advance 
the binary in the galactic potential, and a simple iteration 
scheme to evolve the semi-major ax:is and eccentricity of 
the binary as gravitational radiation drives a and e to zero, 



© 0000 RAS, MNRAS 000, 000-000 



4 Bloom, Sigurdsson, & Pols 



assuming the orbit averaged quadrupole dominated approx- 
imation (Peters 1964) . The integration is continued until ei- 
ther 1.5 X 10^° years have passed (no merger in Hubble time) 
or the characteristic time to merger is short compared to the 
dynamical time scale of the binary in the halo (ie. the binary 
won't move any further before it merges). We then record 
the 3-D position of the binary relative to the presumptive 
parent galaxy and the time since formation. 



4 RESULTS 

4.1 Orbital parameter distribution after the 
second supernova 

Figure ^ shows the distribution of orbital parameters (semi- 
major axis and period) after the second supernova for bound 
NS-NS pairs for the Hansen & Phinney (1997) kick distribu- 
tion (akick = 190 km s~^). As found previously (eg. Portegies 
Zwart & Spreeuw 1996), bound systems tend to follow lines 
of constant merger time. The density of systems in figure 1 
can be taken as the probability density of finding a NS-NS 
binaries directly after the second supernova. In time, the 
shorter-lived systems (higher e and shorter period) merge 
due to gravitational radiation. Thus, at any given time af- 
ter a burst of star-formation there is an observational bias 
towards finding longer-lived systems. In addition, there is 
a large observational bias against finding short period bina- 
ries (Johnston & Kulkarni 1991). That the observed PSR-NS 
systems lie in the region of parameter space with low initial 
probability is explained by these effects. The time-dependent 
probability evolution has been discussed and quantified in 
detail by Portegies Zwart & Yungelson (1998). Figure ^ 
shows the distribution of merger times as a function of sys- 
tem velocity. A majority of systems merge in ~ 10* yr spread 
over system velocities of 50 — 500 km s~^. A subclass of 
systems are have spatial velocity and merger time which are 
anti-correlated. 
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Figure 2. The distribution of merger times after second super- 
novae as a function of system velocity. Left of the vertical line, all 
pairs created are gravitationally bound to a undermassive host 
(3 X 10^^ Mq) at the disk scale radius. Of the pairs that are un- 
bound, only the pairs in the shaded region could travel more than 
~ 25 kpc (linearly) from their birthplace and merge within a Hub- 
ble time (< 1.5 X 10^'^ yrs). Since the spatial velocity of observed 
NS binaries includes both the initial circular velocity of the sys- 
tem and the system velocity due to kicks from each supernova, 
the true system velocities are highly uncertain. For comparison, 
though, we demark the range of accepted kick velocities of PSR 
1913-1-16 with a long rectangle (the merger time is much better 
constrained than that depicted); this illustrates a general agree- 
ment of the system velocity of PSR 1913-1-16 and the modeled dis- 
tribution of bound NS binaries. The slightly longer merger time 
of PSR 1913-1-16 than expected from the density of systems in 
this parameter space is explained in section 4.1 of the text. 



4.2 Coalescence/Birth Rates 

We have explored the consequences of different kick 
strengths on the birthrates of NS-NS binaries. Table 4.2 
summarizes these results. 

Earlier work (eg. Sutantyo 1978; Dewey & Cordes 1987; 
Verbunt, Wijers & Burm 1990; Wijers et al. 1992; Brandt & 
Podsiadlowski 1995) in which asymmetric kicks were incor- 
porated with a single NS component binary (as in LMXBs, 
HMXBs) noted a decrease in birthrate with increased kick 
magnitude. Portegies Zwart & Spreeuw (1996) and Lipunov, 
Postnov, & Prokhorov (1997) found a similar effect on the 
bound NS pair birthrates. Lipunov (1997) provides a good 
review of the expected rates. Clearly the birthrate of NS-NS 
binaries is also sensitive to the total SN type II rate (which 
is observationally constrained to no better than a factor of 
two, and theoretically depends both on the uncertain high 
mass end of the initial mass function and the total star for- 
mation at high redshift), and is also sensitive to the fraction 
of high mass stars in binaries with high mass secondaries. 

We concentrate our discussion of NS-NS binary 
birthrates to galactic systems for which the SN type II is 
fairly well-known (such as in the Galaxy). It is important to 



note, however, that the SN type II rate may be quite high 
in low surface-brightness and dwarf galaxies (eg. Babul & 
Ferguson 1996). This would subsequently lead to a higher 
NS-NS birthrates in such systems than a simple mass scal- 
ing to rates derived for the Galaxy. 

Recently van den Heuvel & Lorimer (1996) find (obser- 
vationally) the birthrate of NS-NS binaries to be 8 Myr~^. 
Lipunov, Postnov, & Prokhorov (1997) find between 100 
and 330 events per Myr in simulations. Portegies Zwart & 
Spreeuw (1996) found birthrates anywhere from 9 to 384 
Myr~^ depending mostly on the choice of asymmetric kick 
strength in their models. We note that our derived birthrate 
of ~ 3 Myr~^ for high (t„ = 270 km s~^ is comparable to 
those found Portegies Zwart & Spreeuw (particularly model 
"ck") with an average 3-D kick velocity of 450 km s~^. Also, 
for low velocity kicks ((Tkick = 95 km s~^) our birthrates ap- 
proach that of Portegies Zwart & Spreeuw models with no 
asymmetric kicks. 

The discrepancies between this and other work, there- 
fore, we believe are largely due to the choices of supernovae 
kick distributions and strengths. That the absolute birthrate 
varies by an order of magnitude depending on the binary 
evolution code and asymmetric kick distributions used in 
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Table 1. The spatial distribution of coalescing neutron stars in various galactic potentials. Though the average distance from center a 
pair travels before coalescence (rfavg) generally decreases with increasing galactic mass, the median distance (limedian) scales with disk 
radius (rdisk) 



Galaxy parameters Coalescence Distance 



Run 


i>circ (km/s) 


»'brcak (kpc) 


''disk (kpc) 


M (lO"Af0) 


L 


(imodian (kpc) 


davg (kpc) 


a 


100 


1 


1 


0.092 


< 0.05L. 


4.3 


66.2 


b 


100 


3 


1 


0.278 


~ O.IL* 


4.0 


50.1 


c 


100 


3 


3 


0.278 


~ O.IL* 


8.7 


68.8 


d 


150 


3 


1 


0.625 


~ 0.5L* 


3.1 


24.8 


e 


150 


3 


3 


0.625 


~ 0.5L* 


7.7 


54.1 


f 


225 


3 


3 


1.41 


~ IL. 


6.0 


29.9 


g 


225 


3 


1 


1.41 


~ IL. 


2.3 


7.1 


h 


225 


5 


3 


2.34 


~ 2Z/^ 


6.0 


21.4 


i 


225 


5 


5 


2.34 


— 27v^ 


9.9 


30.2 



Table 2. The bound NS— NS binary birthrate and merger time properties as a function of supernova kick strength. A Maxwellian 
distribution characterized by a velocity dispersion (<Tkick) is assumed. 



o-kick (km/s) 


Birthrate (Myr ~l) 


Tmcdian (yr) 


T-avg (yr) 


95 


49 


1.4 X 10* 


9.4 X 10** 


190 


10 


7.0 X 10^ 


8.0 X 10* 


270 


3 


5.5 X lO'^ 


7.0 X 10* 



" Average merger time of pairs merging in less than 1.5 X 10^" years. 



different studies, hints at the uncertainty in the knowledge 
of the true birthrates. 



4.3 Spatial Distribution 

Approximately half of the NS-NS binaries merge within ~ 
10* years after the second SN; this merger time is relatively 
quick on the timescale of star-formation. In addition, half 
the pairs coalesce within a few kpc of their birthplace and 
within 10 kpc of the galactic centre (see figure ^ regardless 
of the potential strength of the host galaxy. As shown in 

figure M, galaxies with Mgalaxy > 10^° Mq (L ^ O.IL,), 90 

(95) percent of the NS-NS mergers will occur within 30 (50) 
kpc of the host. In the least massive dwarf galaxies with 
Afgaiaxy ^ 9 X lO^M© (^ O.IL*), 50 (90, 95) percent of 
mergers occur within ~ 10 (100, 300) kpc of the host (see 
figure ^). So, for example, assuming a Hubble constant of 
Ho = 65 km s"^ Mpc"^ and D. = 0.2, we find that 90 (95) 
percent of NS binaries born in dwarf galaxies at redshift z — 
1 will merge within ~ 12.7 arcsec (~ 38.2 arcsec) of the host 
galaxy. These angular offsets can be considered the extreme 
of the expected radial distribution since the potentials are 
weakest and we have not included the effect of projection. 
We would expect 50 (90, 95) percent the mergers near non- 
dwarf galaxies to occur within ~ 1.3 (3.8, 6.4) arcsec from 
their host at 2 = 1 for the cosmology assumed above. 

Given the agreement of our orbital parameter distri- 
bution (figure 1) and velocity distribution (figure 2) with 
that of Portegies Zwart & Yungelson (1998), the discrep- 
ancy between the derived spatial distribution (see figure 8 
of Portegies Zwart & Yungelson) is likely due to our use of a 
galactic potential in the model. This inclusion of a potential 




r (kpc) 

Figure 3. The radial distribution of coalescing neutron stars 
around galaxies of various potentials. The letters refer to runs 
in table 1. In all scenarios, at least 50% of the mergers occur 
within 10 kpc of the host galaxy. The wider radial distribution of 
in the underluminous galaxy scenarios (a,c) reflects the smaller 
gravitational potential of underluminous galaxies. 

naturally keeps merging NSs more concentrated towards the 
galactic centre than without the effect. 
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Figure 4. NS-NS merger rate dependence on redshift. The dot- 
ted line is a reproduction of tlie SFR from Madau 1997 with cor- 
responding units on the left-hand axis. The SFR curve as seen as 
a lower limit to the true star-formation history since dust may ob- 
scure a large fraction SFR regions in galaxies. The right hand axis 
is the (unobscured) GRB rate if the bursts arise from the merger 
of two NS— NS assuming a merger time distribution found in the 
present study (dot— dashed line). Both the SFR and merger rate 
are in co-moving units (assuming Hq = 50 km Mpc~^). The 
normalisation of the burst rate is taken from Wijers et al. 1998. 



5 DISCUSSION 

Although the NS-NS birthrate decreases with increased ve- 
locity kick, the distribution of merger time and system ve- 
locity is not affected strongly by our choice of kick distri- 
butions. Rather, the shapes of the orbital and velocity dis- 
tributions (figures 1 and 2) are closely connected with the 
pre-SN orbital velocity, which is itself connected simply with 
the evolution and masses. That is, bound NS binaries come 
from a range of parameters which give high orbital veloc- 
ities in the pre-second SN system. The orbital parameters 
(and merger time distribution) of binaries which survive the 
second SN are not sensitive to the exact kick-velocity dis- 
tribution. We suspect this may be because bound systems 
can only originate from a parameter space where the kick 
magnitude and orientation are tuned for the pre-second SN 
orbital parameters. The overall fraction of systems that re- 
main bound IS sensitive to the kick distribution insofar as 
the kick distribution determines how many kicks are in the 
appropriate range of parameter space. 

Since NS-NS binaries are formed rapidly (with an aver- 
age time since ZAMS of ~ 22 million years) and the median 
merger time is of order one hundred million years regard- 
less of the kick velocity distribution (see table 2), the rate 
of NS-NS mergers should closely trace the star formation 
rate. In the context of gamma-ray bursts, where merging 
NSs are seen as the canonical production mechanism, this 
result implies that the GRB merger rate should evolve pro- 
portionally to the star formation rate (see figure ^; see also 
Bagot, Portegies Zwart, & Yungelson 1998). Indeed, several 



studies (Totani 1997; Lipunov, Postnov, & Prokhorov 1997; 
Wijers et al. 1998) have consistently fit the GRB log N- 
log P curve to a model which assumes such a rate density 
evolution. 

If indeed gamma-ray bursts arise from the coalescence 
of neutron star binaries, then we confirm that GRBs should 
trace the star formation rate in the Universe; thus most 
GRBs should have redshifts near the peak in star formation 
(currently believed to be 1 < 2 < 2; Madau et al. 1996) al- 
though the observed distribution may be skewed to lower 
redshifts by obscuration at high redshift (eg. Hughes et 
al. 1998). Determination of the distribution of x-ray and 
optical counterparts to GRBs may help constrain the true 
cosmological star formation history, though the observations 
of GRB counterparts are vulnerable to some of the same ex- 
tinction selection effects that complicate determination of 
high redshift star formation rates. Figure ^ illustrates the 
redshift dependence of the GRB rate assuming the bursts 
arise in NS mergers. 

The minimum required local (isotropic) bursting rate of 
0.025 galactic event per Myr (Wijers et al. 1998) is consis- 
tent with our birthrate results (table 2) assuming a beam- 
ing fraction of 1/10-1/100 for the gamma ray emission and 
our canonical values for the type II supernova rate and su- 
pernova binary fraction. The effects of beaming should be 
observed in both the fight curves of GRB afterglow and in 
deep transient searches (eg. Woods & Loeb 1998). 

In the case of GRB 970508, Bloom et al. 1998a and 
Castro-Tirado et al. 1998 (see also Natarajan et al. 1998) 
found that the host is an underluminous dwarf galaxy; the 
close spatial connection (offset < 1") of the OT with the 
galaxy is then a case (albeit weakly) against the NS-NS 
merger hypothesis as the a priori probability is <; 20% (fig- 
ure H) . Paczynski 1998 first pointed out that the close spa- 
tial association with a dwarf galaxy is a case against the 
NS-NS merger hypothesis. Certainly more transients are 
required to rule against the NS-NS merger hypothesis; we 
note, however, that dust obscuration and projection effects 
may severely bias the sample (see above discussion). 

The verdict on the reconciliation of the expected radial 
distribution of NS-NS mergers with hosts of known GRBs is 
still out. Sahu et al. 1997 found the optical transient associ- 
ated with GRB 970228 to be slightly offset from the centre 
of a dim galaxy but without a redshift it is still unclear 
as to the the true luminosity of the host and thus the ex- 
pected offset of the OT in the NS-NS merger hypothesis. 
Similarly, small or negligible offsets of GRB afterglow with 
faint galaxies has been found for GRB 980326 and GRB 
980329 (Djorgovski et al. 1998). Kulkarni et al. 1998 found 
the redshift of the purported host galaxy of GRB 980329 to 
be z = 3.4 implying the host is L J> L»; the expected offsets 
of NS-NS mergers around massive galaxies (figure ^, models 
d through i) is then consistent with their finding of an OT 
offset ~ 0.5 kpc. 

A few well-established offsets cannot tell us what is the 
true distribution of GRBs around host galaxies. As more 
OTs are discovered, we will hopefully build up a large sam- 
ple to statistically test statistically the offsets. Fortunately 
the unobscured afterglow emission strength is coupled with 
the density n of the surrounding interstellar medium (ISM) 
with intensity scaling as ^Jn (Begelman, Meszaros & Rees 
1993; Meszaros, Rees, & Wijers 1998); however, high ISM 
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densities tracing dust will tend to obscure rest-frame UV 
and optical emission from the transient. In the absence of 
strong absorption from the surrounding medium transients 
of GRBs are preferential found close to where they are born, 
in the disk. However, dust obscuration and projection ef- 
fects severely complicate determination of the true offset 
of OTs from their host galaxy. Furthermore, identification 
of the host with a GRB becomes increasingly difficult with 
distances beyond a few light radii (~ 10 kpc) of galaxies 
(although see Bloom, Tanvir & Wijers 1998). 

If all afterglows, especially those where little to no ab- 
sorption is implied, are found more highly concentrated than 
predicted in figure ^, the NS-NS merger hypothesis would 
lose favour to models which keep progenitors more central 
to their host. GRBs as events associated with single massive 
stars such as microquasars (Paczynski 1998) or failed type 
lb SN (Woosley 1993) could be possible. Alternatively, one 
may consider neutron star-black hole (BH) binaries as the 
progenitors of GRBs (Mochkovitch et al. 1993; Meszaros & 
Rees 1997). Most black hole X-ray binaries have low-spatial 
velocities (although Nova Sco has Vsys — 100 km s~^; see 
Brandt, Podsiadlowski, & Sigurdsson 1995) so NS-BH bi- 
naries should have system velocities ~ 3 to 10 times smaller 
than NS-NS binaries. One would expect NS-BH systems 
to be borne with higher eccentricities than NS-NS systems 
leading to quicker merger. Moreover, NS-BH binaries are 
more massive than NS-NS binaries and merger time due to 
gravitational radiation scales strongly with mass. Thus the 
attraction is that NS-BH mergers would be preferentially 
closer to their host and their merger rate might be small 
enough so as to require no beaming. Alternatively, gamma- 
ray bursts could arise from several of these plausible pro- 
genitor models and still be consistent with basic relativistic 
fireball models. 



6 CONCLUSIONS 

A reconciliation with the expected distribution of presumed 
progenitors of GRBs and observed transient/host offsets is 
clearly required. We find for all plausible galactic potentials 
that the median radial offset of a NS-NS merger is less than 
10 kpc. And in all but the most shallow potentials, ninety 
percent of NS-NS binaries merge within 30 kpc of the host. 
At a redshift of z = 1 (with Ho — 65 km s~^ Mpc~^ and 
— 0.2), this means that ninety percent the coalescences 
should occur within ~ 4 arcsec from the host galaxy. Al- 
though the expected spatial distribution of coalescing neu- 
tron star binaries found herein is consistent with the close 
spatial association of known optical afterglows of gamma-ray 
bursts with faint galaxies, a non-negligible fraction (~ 15%) 
of GRBs should occur well outside dwarf galaxy hosts if the 
NS-NS hypothesis is correct. Otherwise, other models which 
keep progenitors closer to their host (eg. BH-NS mergers, 
microquasars, or "failed SN type lb") would be preferred. 

As all the progenitor models mentioned are connected 
with high-mass stars, the true GRB afterglow rate as a func- 
tion of redshift should trace the star-formation rate in the 
Universe. However, environmental effects, such as dust ob- 
scuration, may severely bias the estimate of the true offset 
distribution. Even in the NS-NS models where progenitors 
have a natural mechanism to achieve high spatial veloci- 



ties, most will be closely connected spatially to their host. 
Redshifts derived from absorption in the afterglow spectra 
should be nearly always that of the nearest galaxy (Bloom 
et al. 1998b). Rapid burst follow-up (< 1 hr), with spectra 
taken while the optical transients are bright should confirm 
some form of absorption from the host galaxy. 

We have confirmed the strong dependence of birthrate 
of NS-NS binaries on kick velocity distribution and found 
the independence of the orbital parameters after the second 
supernova (and hence merger times and spatial velocity) on 
the choice of kicks. The methodology herein can be extended 
to include formation scenarios of black holes. This could pro- 
vide improved merger rate estimates for LIGO sources, and 
estimate the relative contribution of coalescences between 
neutron stars and low mass black holes to the event rate. 
Detailed modeling of the Milky Way potential would also 
allow predictions for the distribution of NS-PSR binaries 
observable in the Milky Way, which would provide an inde- 
pendent test of the assumptions made in these models. 
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